Generation of strong magnetic fields by r-modes in milfisecond accreting neutron 
stars: induced deformations and gravitational wave emission 
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DifTerential rotation induced by tlie r-mode instability can generate very strong toroidal fields 
in the core of accreting, millisecond spinning neutron stars. We introduce explicitly the magnetic 
damping term in the evolution equations of the r-modes and solve them numerically in the Newtonian 
limit, to follow the development and growth of the internal magnetic field. We show that the 
strength of the latter can reach large values, B ~ 10^* G, in the core of the fastest accreting 
neutron stars. This is strong enough to induce a significant quadrupole moment of the neutron star 
mass distribution, corresponding to an ellipticity |eB| ~ 10~*. If the symmetry axis of the induced 
magnetic field is not aligned with the spin axis, the neutron star radiates gravitational waves. We 
suggest that this mechanism may explain the upper limit of the spin frequencies observed in accreting 
neutron stars in Low Mass X-Ray Binaries. We discuss the relevance of our results for the search 
of gravitational waves. 

PACS numbers: 04.40.Dg, 04.30.Tv, 04.30.Db, 97.60.Jd 



I. INTRODUCTION 

The r-mode instability plays an important role in the 
physics of millisecond neutron stars NSs. It excites the 
emission of gravitational waves (GWs), which carry away 
spin angular momentum causing the star to spin down. 
It also gives rise to large scale mass drifts, particularly in 
the azimuthal direction, and to differential rotation [1-5] . 
Differential rotation in turn produces very strong toroidal 
magnetic fields inside the star and these fields damp the 
instability by tapping the energy of the modes. This 
mechanism has been investigated in the case of rapidly 
rotating, isolated, newly born neutron stars in Refs. [1- 
3, 6] and in the case of accreting millisecond neutron and 
quark stars in Refs. [7, 8]. Magnetic fields deform the 
star and if the magnetic axis is not aligned with the ro- 
tation axis the NS undergoes free body precession. The 
deformation induced by the strong toroidal field is such 
that the symmetry axis of the precessing NS drifts on a 
timescale determined by its internal viscosity, eventually 
becoming an orthogonal rotator [9]. This is an optimal 
configuration for efficient GW emission, thus enhancing 
angular momentum losses from the NS. 
In this work we resume the analysis in Ref. [7] about the 
growth of the internal magnetic fields in accreting mil- 
lisecond neutron stars induced by r-modes. We account 
for the back-reaction of the generated magnetic field on 
the mode amplitude, by introducing the magnetic damp- 
ing rate into the evolution equations. 
We find that the GW emission due to the magnetic de- 
formation grows secularly in accreting NSs, effectively 
limiting the growth of their spin frequency. In particu- 
lar, this mechanism could naturally explain the observed 



cut-off above 730 Hz in the distribution of spin frequen- 
cies of the faster accreting neutron stars in Low Mass 
X-Ray Binaries (LMXBs) [10, 11]. 

Finally, we estimate the strain amplitude of the GW sig- 
nal emitted from LMXBs in the local group of galaxies 
to establish their detectability with advanced LIGO and 
Virgo interferometers. 

The paper is organized as follows: in Sec. II we write 
the r-mode equations taking into account also the mag- 
netic damping rate; in Sec. Ill we discuss the magnetic 
deformation of the star and the associated gravitational- 
waves emission; in Sec IV we derive the expression of the 
magnetic damping rate; in Sec. V we discuss numeri- 
cal solutions of the equations; in Sec. VI we deal with 
the evolutionary scenarios for accreting neutron stars; in 
Sec. VII we focus on the relevance of our results for the 
search of gravitational waves. Finally in Sec. VIII we 
summarize our conclusions. 



II. R-MODE EQUATIONS 

We derive here the equations that describe the evolu- 
tion of r-modes in a NS core with a pre-existing poloidal 
magnetic field, Bp, and the associated generation of an 
internal azimuthal field B^. Following Ref. [12], the to- 
tal angular momentum J of a star can be decomposed 
into an equilibrium angular momentum J^ and a canon- 
ical angular momentum Jc proportional to the r-mode 
perturbation: 



J ^J^{M,n) + {l-Kj)J^ 



Jc = -K^a^J, (1) 
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where Kij^) ^re dimensionless constants and J* = 7*0. 
The canonical angular momentum obeys the following 



equation [7, 13]: 

dJc/dt^ 2Jc{F;^{M,n) 



iF,iM,n,n) + F^^{M,n,B)]} (2) 



where F, 
the 



is the gravitational radiation growth rate of 
r-mode, F^ = Fs + Ff, is the sum of the shear and 
bulk viscous damping rate and Fmi is the damping rate 
associated with the generation of an internal magnetic 
field. Finally, Ty{t) is a spatially averaged temperature. 
The total angular momentum satisfies the equation: 



dJ/dt ^ 2JcFl + Ja{t) - FnF„ 



(3) 



where Ja is the rate of variation of angular momentum 
due to mass accretion (we assume Ja = M{GMRy^'^, see 
Ref. [14]) and Fm^ is the magnetic braking rate associ- 
ated to the external poloidal magnetic field. In this work 
we do not consider additional torques such as, e.g. the 
interaction between the magnetic field and the accretion 
disk. Combining Eqs. (2) and (3) we obtain the evolution 
equations of the r-mode amplitude a and of the angular 
velocity of the star fl: 



da 



a{Fl-F,-F^,) 

+a[K,Fj; + (1 - K,){F, + Fm,)]K, 



aM 



G 



aF„ 



dn 

'dt 



2m \MR3J ' 2 
-2K,na^[KjF;^ + (1 - Kj){F, + F„J] 

"^T \MR3 



-nF„ 



(4) 



(5) 



with / — 0.261 for an n=l polytrope 
^^, see Ref. [15]. Our results turn out 



M 

where /, = IMR'^ 
and Kc ^ 9.4 x 10" 
to be rather insensitive to the exact value of Kj ^ 1 (see 
Ref. [12]). 

The classical r-mode equations are simply recovered from 
the above by neglecting the magnetic damping term, i^„i. . 
In this case, it is readily seen that the r-modes will grow 
only if the condition F^ > F^ is satisfied, while viscosity 
keeps the NS stable against r-modes if this condition is 
reversed. This simple argument defines the classical in- 
stability window for r-modes, which is drawn in Fig. 1 as 
the thick dot-dashed curve. 

Well inside the classical instability window, where vis- 
cosity is negligible, a new type of equilibrium can be 
reached. Eq. (4) and Eq. (5) indeed imply that the 
growth of r-modes can be quenched by the magnetic 
damping term, F^^, when the condition F^ < F^m holds. 
Since F^ — F^{fl) and i^m. = Fjn.{fl, B), the equilibrium 
condition FJ" = F^m implies a relation il = ^{B), which 
can be written as: 



68 Bli' B'{1 Mi.4 



2/7 n-3/7 



R 



10 



(6) 



where Bd.s = 

(M/1.4Mo) 



: (Bd/IO^G), B^,i2 = (B^/lQi^G), Mi.4 
R^Q = (i?/10km). 



Here, and in the following, i/ indicates the NS spin fre- 
quency and i?0 is an energy averaged, azimuthal mag- 
netic field induced by the r-modes. 

For a given combination of (Bd, B^), Eq. (6) represents 
the critical frequency above which the star is r-mode un- 
stable even against the effect of the magnetic damping 
term. For a given value of the spin frequency, on the 
other hand, Eq. (6) tells us the minimal combination of 
field strengths required to damp the r-modes. With the 
value of Bd fixed, r-modes are excited as long as B^ is 
lower than required by Eq. 6. They continue to generate 
further azimuthal field, until Eq. (6) is satisfied and the 
NS is eventually stabilized with respect to r-modes. 



III. MAGNETIC DEFORMATION AND GW 
EMISSION 

We assume a NS with an internal poloidal field which 
is described by Ferraro's solution [7, 16] 



B" 
Bd 



'(t = 0) 



(7) 
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where Bd is the surface field strength at the stellar equa- 
tor. This is matched to an exterior dipolar field aligned 
with the star spin axis 



B" 



i?3 
= Bd-5-(2cos( 



smt'ee) 



(8) 



As for the internal structure of the NS, it is widely ac- 
cepted that protons will be superconducting at tempera- 
tures T < 10^ K [17, 18], at least in a fraction of the NS 
core. To maintain focus on the salient properties of our 
scenario, we carry out here a detailed study of the growth 
of r-modes and of the azimuthal field in a normal fluid 
core. In Section VI B we consider the likely occurrence of 
a shell of superconducting protons in the outer core and 
assess its impact on our conclusions. 
Once modes are excited at an initial time to i they induce 
azimuthal drift motions of fluid parcels in the NS core. 
Following Ref. [2], the total angular displacement from 
the onset of the instability up to time t reads 



A#(r,t)-^(-^)fc2(0)y a\t')n{t')dt' + 0{c 



(9) 



where k2{0) = (l/2)^(5!/7r)(sin26l - 2cos26l). Magnetic 
held lines are twisted and stretched by the shearing mo- 
tions and a new held component, in the azimuthal direc- 
tion, is generated accordingly. The relation between the 
new and the original magnetic flcld inside the star in the 
Lagrangian approach is [2] : 



B^ B^ dx it] 

— (x,t)^— (x,to) ^ l;' 



(10) 
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FIG. 1: R-mode instability regions for NSs with strong in- 
ternal fields. The r-mode instability window without internal 
magnetic field corresponds to the thick dotted-dashed line. 
The vertical dashed lines indicate the internal equilibrium 
temperatures of accreting neutron stars in LMXBs for, re- 
spectively, a maximum mass accretion rate M = f O~*M0yr~^ 
(right) and a minimum value of M = lO~^^M0yr~^ (left). 
The instability window with a stable configuration of mixed 
poloidal-toroidal fields in the inner core, obtained from the 
condition Fg — Fm^ + F,^, is indicated by thin dotted-dashed 
lines. Three different values of the internal magnetic field are 
considered and we assume that poloidal and toroidal compo- 
nents have similar strengths. The light blue area indicates the 
region accessible to stable accreting NSs (see Sec. VI). The 
thick solid line indicates the maximum spin frequencies when 
the emission of gravitational waves due to the star's deforma- 
tion is taken into account {B^ ~ 10^* G, |es| = 10~*). 



This equation implies that the radial dependence of the 
initial and final magnetic field is the same. Integrat- 
ing over time the induction equation in the Eulerian ap- 
proach one gets [3]: 

SB'^ ~ 55'' ~0 (11) 

5B'^{r,t) ~ B^ f <j){t')dt' ~B^oAS:'''{r,t) (12) 

where B"^ is the toroidal component generated by the r- 
modes. 

A toroidal field B^ deforms a neutron star into a prolate 
ellipsoid with ellipticity 



/. 



h 



£B 



(13) 



where Ijk = Jy p{r){r'^Sjk — XjXk)dV is the inertia tensor. 
For a neutron star with a normal core, this gives 



es = -h X 10^12 ^4^ jy,j~2 



4 -°<*,12 



(14) 



where fee ~ [0.5-5] is a parameter depending on the con- 
figuration of the internal fields [9, 19, 20]. 



In general, we expect the symmetry axis of the magnetic 
deformation not to be perfectly aligned with the rotation 
axis. In this situation, J* precesses around the symmetry 
axis with a period Pproc — P/\£b\ [21^23], P being the 
spin period of the star. 

Due to the neutron star rotation, the relaxed state of 
its crust also has an intrinsic oblate deformation. This 
corresponds to the shape it would retain if the NS spin 
were stopped completely, without the crust cracking or 
re-adjusting in anyway. The ellipticity associated to this 
crustal reference deformation can be written as [9, 24, 25]: 



'"">-' irni)' 



(15) 



where kc — [1-3] for an accreted crust [25]. As long as 
this is dominant, the NS effective mass distribution will 
be that of an oblate ellipsoid. In this case, the spin and 
symmetry axis of the freely precessing NS will tend to 
align. However, as the internal toroidal field grows, the 
ensuing magnetic deformation will at some point over- 
come the maximal oblateness that the crust can sustain 
(Eq. (15)). The NS first enters the r-mode instability 
window when its spin rate is t's,! ^ 150 Hz. Therefore, the 
magnetic deformation will first become dominant when 



leBJ > 2x10 



-10 



J^s.i 



150 Hz 



es. 



B^ > 1.5xlO^^G( ^ 



1/2 



150 Hz 



Mi.4 i?ro' -(16) 



Beyond this point, the total NS ellipticity becomes domi- 
nated by the magnetically-induced, prolate deformation. 
The freely precessing NS now becomes secularly unsta- 
ble. In the presence of a finite viscosity of its interior, 
the wobble angle between the angular momentum J^, and 
the symmetry (magnetic) axis grows until the two be- 
come orthogonal. This occurs on a dissipation (viscous) 
timescale, Tv, which can be expressed as [9] 



nPn 



10 



11 



1 105 J 



100 Hz 



10- 



|es| 



(17) 

The parameter n measures the dissipation timescale in 
units of the free precession cycle. Below, we determine 
quantitatively its role in promoting the orthogonaliza- 
tion of the NS symmetry and rotation axes. Given the 
significant theoretical uncertainties on the value of n, we 
normalize it to an educated guess [9] but will consider it 
effectively as a free parameter. 

An orthogonal rotator with a time- varying quadrupole 
moment, Q22 ~ I*^b^ emits gravitational- waves at a 
rate [26] 



S^ = - 

gw 



32G{heBf 



n' 



(18) 



which produces spin down at the rate _F? = —{E^^/2E) 
1 - „ / f^o \ 2 



» 5 X 1014 



h Ml. 4 R 



Vio-sy 



p_t s-i 



(19) 



where h = //0.261. 

Therefore, if the orthogonahty condition is satisfied, the 
NS looses significant spin angular momentum to GWs 
also through the induced magnetic deformation. This is 
accounted for by introducing the appropriate GW radia- 
tion term, _F?, in the Eq. (3), which now becomes: 



dJ/dt = 2JcF; + Ja{t) - IMFJ^ + F, 



(20) 



Eq. (4) and Eq. (5) for the r-mode amplitude and the NS 
spin frequency are modified accordingly: 



of a in detail, we note that a new asymptotic equilib- 
rium could take place in this situation, with the material 
torque being balanced by the GW torque due to the mag- 
netic deformation, F^ . This new equilibrium condition, 

J a = t/(f , leads to a limit spin frequency 
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t^ma. «990/i ''^Mi.4'"'i?10 
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Eq. (21) and Eq. (22) now contain two terms for GW 
emission. Of these, F^ is given by Eq. (19) while we 
adopt for F^ the r-mode gravitational radiation reaction 
rate due to the I = m = 2 current multipole [27] 



IV. MAGNETIC DAMPING RATE 

Following Ref. [3], we make use of Eq. (9) and Eq. (12) 
to derive the expression for the magnetic damping term. 
Let us first write the variation of the magnetic energy for 
a NS with a normal fluid core 



5Em = ^ 



SB^dV 



(27) 



where 6B^ = (SBpf + (5B^f ~ {5B^f = Bj. In the 

last steps wc have taken into account the Eqs. (11,12). 
The rate of variation of magnetic energy can be obtained 
from the Eqs. (9,12) and reads 



^; = -, M,4 



R'^ 



p-6 



(23) 



The evolution equations (21,22) will hold only after 
the symmetry axis of the azimuthal field has become or- 
thogonal to the spin axis. 

If t measures the time since the generated magnetic field 
exceeds the critical value given by Eq. (16), the rotation 
of the magnetic axis will occur when t > Tv(i), where 
Tv is given by Eq. (17). Once the condition t > Tv(t) is 
met and Eq. (21) becomes effective, it is relevant to know 
which of the two GW emission terms is dominant. To see 
this, we consider the ratio between the GW torque due 
to r-modes, Jg — 2JcFg, and the GW torque due to the 
magnetic deformation J^ = I^F^ (see Eq. (20)) 
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VlO-8y' 
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(24) 



The magnetic deformation becomes dominant only when 
the r-mode amplitude gets very small, and the internal 
magnetic field accordingly large. In particular, the ra- 
tio in Eq. (24) becomes smaller than unity once a has 
decreased below the critical value: 



< 2 X 10"' i? 
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(25) 



A numerical solution of Eqs. (21,22) is required to verify 
when the condition given by Eq. (25) is satisfied during 
the system's evolution. Before calculating the evolution 



SE 



2ABJR^ 
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x2 

a\t')n{t')dt' . (28) 



From this the expression of the magnetic damping rate 
is derived [7]: 



^^ ^ (dE^/dt) 
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Mn 



(29) 



where E ~ 8.2 x 10^ '^a'^Mfl'^R'^ is the energy of the r- 
mode [3] and A ~ 0.99 [7]. We make use of Eqs. (27) and 
(28) to estimate the toroidal magnetic field generated by 
r-modes 



B, 



4A 
3^ 



1/2 



Bd I a^{t')n{t')dt' . (30) 
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Note that F^j is not related to a in a simple way. The 
magnetic damping rate depends, like Bj, on the inte- 
grated history of the r-mode amplitude since the star 
first entered the instability window (cfr. Eq. (12)). 
Finally, by using the volume-averaged expression for the 
azimuthal field, B^, we can re- write the expression of the 
magnetic damping rate as 



F 



'N 



^^^—^^R,,M^lBa^,B^ 
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F-3S- 



(31) 



V. NUMERICAL SOLUTIONS 

In the scenario we are describing, an initially slowly ro- 
tating NS is secularly spun up by mass accretion. When 
its spin frequency reaches a few hundred Hz, the NS en- 
ters the classical r-modes instability window. As the in- 
stability develops, the evolution of the r-mode amplitude 
becomes coupled to the growth of an internal, toroidal 
magnetic field. We now turn to a detailed numerical cal- 
culation of their coupled evolution, in the light of our 
previous discussion. This will help us clarify the se- 
quence of events expected to occur as an accreting NS 
in a LMXB is secularly spun up by the material torque. 
In all our calculations we consider values of the mass 
accretion rate M = (10^^-10"^°) Mgyr'^ Note that, 
in accreting NSs, the mass accretion rate has an upper 
limit MEdd ^ lO^^Mgyr^i and most LMXBs do not ac- 
crete at this rate for a long time. Here M-^dd is the mass 
accretion rate that produces the Eddington luminosity. 



A. Numerical estimates of relevant rates 

We begin by discussing the main physical quantities in 
the evolution equations that were not described in previ- 
ous sections. 

For non-superfluid matter the shear viscosity damping 
rate reads [271 



F.,= 



1 



6.7 X 10^ 



,,^5/4 p-23/4„-2 ^- 



(32) 



where Tg = T/10^ K, while the bulk viscosity damping 
rate is given by [15] 



Fh^ 
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which we can approximately rewrite as 



(33) 



(34) 



However, for the temperatures of interest here the bulk 
viscosity damping rate is a few orders of magnitude 
smaller than the shear viscosity damping rate and there- 
fore it is negligible. 

Notice that viscous damping of the modes depends 
strongly on temperature, and that temperature will in 
turn be affected by viscous heating. It is thus impor- 
tant to describe accurately the global thermal balance of 
the NS. We consider three main factors: modified URCA 
cooling (cu), shear viscosity reheating (eg) and accretion 
heating (e„). The equation of thermal balance of the star 
therefore reads: 



d 
dt 



1 



-C,T 



^u ~r ^s I ^T) 



(35) 



Here, Cv is the total heat capacity of the NS [28]: 

Cv = 1.6 X lO^^Mi/^Tg erg K'K (36) 



The cooling rate due to the modified URCA reactions, 
iu, reads [26] 



7.5 X lO^^MfZ/Tl 



erg s 



(37) 



If direct URCA processes can also take place, the NS 
cooling rate is expected to become much larger [29] . A re- 
cent phcnomenological analysis of their implications can 
be found e.g. in Ref. [30]. We assume, for simplicity, that 
only modified Urea processes are allowed in the NS core. 
Dissipation of the r-mode oscillations by the action of 
shear viscosity will contribute to heating the NS. The 
corresponding heating rate, eg, reads [27] 



2o?Vl?MR^JFs 



= 8.3 X lO^^a^f]^ j^9/^4^-^i5/4y-2 ^^g g-i (33) 

where J = 1.635 x 10"^^ 

The accreting material exerts an extra pressure onto the 
NS crust causing direct compressional heating. The com- 
pression also triggers pycnonuclear reactions in the crust, 
which release further heat locally. The total heating rate, 
given by the sum of the two contributions, is [31] : 

£„== X 1.5 MeV = 4x 10^^Mi4 erg s"^ (39) 

niN 

where tojv is the mass of a nucleon and Mi. 4 — M/IAMq 
is measured in s~^. 

We can now solve self-consistently Eqs. (4,5,29,35), to 
obtain the evolution of the core temperature, r-mode am- 
plitude and generated internal magnetic field, -B^. We 
show in Fig. 2 the evolution of the r-mode amplitude, a, 
after the star first enters the classical r-mode instability 
window. Three different values of accretion rate M and 
two values of the initial poloidal magnetic field Bd are 
considered. We find that the maximum values of a are 
in the range amax ^ [10^® - 10"'']. 

In Fig. 3 we show the temporal evolution of the gen- 
erated, volume averaged toroidal magnetic field. Also 
plotted are lines along which the condition Tv(t) — t is 
met, for different values of the parameter n. The secular 
velocity field associated to the r-modes, which is predom- 
inantly azimuthal, clearly induces very large secular ef- 
fects. In particular, very strong toroidal magnetic fields, 
B'^ ^ [10^^-10^''] G, can be produced by the wrapping 
of the pre-existing poloidal field lines, in 10'' — 10^ yrs. 

The implications of these results for the spin equilib- 
rium condition (Eq. (26)) will be discussed in the next 
section. 



VI. THE OVERALL EVOLUTIONARY 
SEQUENCE: A GENERAL DISCUSSION 

In the following we describe likely evolutionary scenar- 
ios for accreting ms spinning NSs, based on the numerical 
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FIG. 2: Temporal evolution of the r-mode amplitude 
a leaving out the instabilities of the magnetic fields. 
We consider three values of mass accretion rate M — 
(lO"*,lO"^,lO"^°)M0yr"^ and two values of the initial 
poloidal magnetic field: Bd ~ 10* G (upper panel) and 
Bd = ICP G (lower panel). 



solutions and the physical discussion of the previous sec- 
tions. We also consider implications of the possible oc- 
currence of the Tayler instability in the strongly twisted 
internal field. Finally, the occurrence of proton super- 
conductivity in (at least a fraction of) the NS core and 
its impact on our arguments are discussed. 

Note that our analysis does not apply to compact stars 
in which exotic matter, like e.g. hyperons or deconfined 
quarks, is present in the inner core. In this type of stars 
the viscosity and thus the r-mode instability window is 
quite different from that of normal NSs [32]. Moreover, 
the evolution of the internal magnetic fields and the mag- 
netic properties of the compact star are also expected to 
be different [8, 33]. 



A. Normal Core 

The growth rate of the azimuthal magnetic field once 
the accreting NS enters the r-mode instability window 
depends on the strength of the initial poloidal field Bd, 
and on the mass accretion rate M . 

In the most favorable cases, e.g. for Bd.% = 1 and 
Af_8 — (0.01 — 0.1), the magnetic field can reach a 
huge strength, B^ « 10^"* G, in a few hundred years. 
In other cases, the internal magnetic field evolves on a 
much longer timescale up to million years (see discussion 
in Ref. [7]). 

As already stated in Sec. Ill, the NS enters the r-mode 
instability window at i^s;! ^ 150. Hence, a minimal mag- 
netic distortion |eB,min| ~ 2 x 10~^° must be achieved 
for the magnetic deformation to overcome the crustal 
oblateness. This corresponds to a minimal azimuthal 



FIG. 3: Temporal evolution of the generated toroidal 
magnetic field obtained by solving Eqs. (4,5). We 

considered three values of mass accretion rate M = 
(10"*, 10"^, 10"^°) Mq yr"^ and two values of initial poloidal 
magnetic field Bd = (10*, 10^) G. The black dashed fine indi- 
cates the minimal strength of the generated toroidal magnetic 
field for which the magnetic deformation becomes dominant 
(Eq. (16)). The time at which the rotation of the internal field 
takes place (i = Tv(i)) is shown as a function of time for five 
values of the parameter n = [10*-10*] (dashed-dotted gray 
lines). Except in cases with large values of initial poloidal 
fields Bd ~ 10^ G and small values of mass accretion rate 



M < 10"® M0 yr"\ huge magnetic fields B"^ 
are generated in less than 10^ years. 



[10"-10'^] G 



field B, 



0,inin 



10 G. Only once this field strength is 



reached the star becomes secularly unstable. However, 
the magnetic axis is driven orthogonal to the spin axis in 
a time Tv (Eq. (17)). Therefore, the condition Tv(i) < t 
must also be satisfied for tilting of the magnetic axis to 
be effective. At the time at which both requirements are 
met, call it tort, the toroidal magnetic field will have an 
intensity i?0,ort and, of course, rv(iort) < ^ort- 

After tilting, the original azimuthal field i?0,ort will 
have acquired an r ~ 9 (poloidal) component, whose 
strength will be comparable to that of the new (f>- 
component. We assume that i3p,ncw = -B^.now ~ 
{l/^/2)Bcf,^oi:t, although the exact ratio will depend on 
the details of the distribution of the internal field. 
Note that the generated magnetic fields may directly af- 
fect the r-mode oscillations when B ^ lO""^^ G [6, 34]. 
We do not include this complicated back-reaction in our 
calculations. Our assumption will be justified a posteri- 
ori, since the maximum generated magnetic fields in our 
model reach values of order 10^"* G. 
In the following we will outline two possible evolution- 
ary paths: which of the two will be followed by the star 
depends on the different combination of the parameters 
(n,fcg,fcc), on the mass accretion rate M and of the initial 
magnetic field Bd- 



1. Evolutionary Path (B^^ort > B^^min) 

The first evolution path we discuss takes place if the 
condition i > Tv is met when B^oit ^ -B^.min- In this 
case a very large deformation of the NS can be obtained. 
The main dissipation mechanism of these magnetic fields 
in the normal fluid core of neutron stars should be am- 
bipolar diffusion whose timescale is given by [35] 



t. 



amh ■ 



3x 10 



rp2T2 



yr 



(40) 



where L is the size of the region embedding the magnetic 
field and L^ = L/IO^ cm. Note that if B^.ort > 10" G 
the generated internal field would decay on a time shorter 
than a few million years to strengths of the order of 
10^^ G which correspond to a deformation jesl ^ 10^^. 
The star is now stabilized with respect to r-modes up to 
frequencies v^r ~ 950 Hz (see Eq. (6)), but the maximum 
spin frequency is limited by GW emission due to the mag- 
netic deformation (see Eq. (26)). The limiting frequency, 
for typical values of M ^ lO^^Afoj/r^^, corresponds to 



630 lil^J^ {€b/10' 



'2/5 Hz. 



2. Evolutionary Path (B^.ort ~ B^^min) 



The second evolutionary path takes place if the condi- 
tion t > Tv is reached when Bcj)^oit ~ B^^^in- After the 
rotation of the magnetic axis, the GW emission term F?' 
becomes effective, with |eB| ~ eB,mm. The new magnetic 
field will have poloidal and toroidal components of com- 
parable strength which stabilize the star against r-modes 
up to frequencies I'cr ~ 500 Hz. Up to this frequency the 
magnetic deformation will thus be the only effective cause 
of GW emission. Mass accretion will continue to spin up 
the star, which may eventually enter again the r-mode 
instability window thus starting the generation of new 
azimuthal field. In Fig. 4 we show the temporal evolution 
of the generated toroidal field and of the r-mode ampli- 
tude a (obtained by solving Eqs. (21,22)) when the star 
enters again the instability window at Vg ~ 500 Hz. The 
magnetic fields can grow further and the star may be sub- 
ject to a new instability when the magnetic field exceeds 
a value of the order of 10^"* G at frequencies I's ~ 650 Hz 
(see Eq. (16) and Fig. 4). After the second flip the evolu- 
tion becomes similar to that described in the first scenario 
being again B^^ort ~ 10^** G and jesl « 10^®. 



Tayler Instability 

The issue of the stability of the magnetic field gener- 
ated by r-modes still needs to be addressed. In the sta- 
bly stratified environment of a stellar interior the Tayler 
instability (or pinch-type) is driven by the energy of the 
toroidal magnetic field. In a slowly rotating NS the insta- 
bility is expected to set in at field strength B^^ci- ^ 10^^ Q 
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FIG. 4: Upper panel: Temporal evolution of the generated 
toroidal magnetic field when the star enters again the r-mode 
instability window at frequencies i/cr ~ 500 Hz after the flip 



of the internal magnetic field {B, 



(^.ort 



BrI, 



10" G). We 



show also the spin frequencies of the star (dashed black lines) 
obtained by Eq. (6). The red dashed line indicates the mini- 
mal strength of the generated toroidal magnetic field to have 
a new instability by magnetic deformation (Eq. (16)). Bottom 
panel: Temporal evolution of the r-mode amplitude a after 
the flip of the internal magnetic field. We consider three val- 
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FIG. 5: Upper panel: Temporal evolution of the generated 
toroidal magnetic field after the development of the Tayler 
instability. Here we assume Bd — 10^^ G. We show also the 
spin frequencies of the star (dashed-dotted lines) obtained by 
Eq. (6) . The red dashed line indicates the minimal strength of 
the generated toroidal magnetic field to have a new instability 
by magnetic deformation (Eq. (16)). Bottom panel: Temporal 
evolution of the r-mode amplitude a after the development 
of the Tayler instability. We consider three values of mass 
accretion rate M = (10"*, lO"'', 10"^°) Mq yr"\ 



[7, 36], with unstable modes growing approximately on 
the Alfven time-scale as long as the ratio of the magnetic 
energy Eem to rotational kinetic energy Trot 
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5 ms 



(41) 



is greater than 0.2 [37]. 

In systems of interest to us, however, this ratio is ex- 
pected to be extremely small, as shown by the estimate 
of Eq. (41). In this condition it is not clear at all that 
the Tayler instability will ever set in [37-39] . It is inter- 
esting to briefly sketch the implications of this instability 
for the evolution of the magnetic fields in fast accreting 
compact stars. 

The most important property of the Tayler instability in 
this context is that, after it develops, the toroidal compo- 
nent of the field produces, as a result of its decay, a new 
poloidal component which can itself be wound up, clos- 
ing the dynamo loop [40] . When the differential rotation 
stops, the field can evolve into a stable configuration of 
a mixed poloidal-toroidal twisted-torus shape, with the 
two components having a comparable strength [41-44]. 
Here we assume that the Tayler instability develops in 
the Alfven time-scale when B^ ^ 10^^ G. At that point a 
poloidal field is generated, having a comparable strength, 
and the star is stable against r-modes only for frequen- 
cies slightly larger than 200 Hz. Mass accretion acceler- 
ates the star back into the instability region and a new 
toroidal field can then develop. The growth of the modes, 
and the generation of the new toroidal field, will start at 
a much higher initial poloidal field, say Bd ~ 10^^ G. 
In Fig. 5 we show the evolution of the r-mode amplitude 
and of the magnetic field in this scenario, after the star 
re-enters the instability window. 

The magnetic fields can grow further and the star may 
be subject to magnetic instability when B^f, > 10^^ G at 
frequencies i^g ~ 450 Hz (see Eq. (16) and Fig. 5). Af- 
ter the flip of the internal magnetic field the evolution 
becomes once again similar to that described in the first 



scenario with _B, 



0,ort 



10" G and [esl « 10" 



B. Superconducting layer 

Protons in a NS core are expected to undergo a transi- 
tion to a (type-II) superconducting state at temperatures 
< lO'' K [45]. Recent calculations of the proton energy 
gap indicate that this transition should occur only in a 
limited density range in the outer core [17, 46], corre- 
sponding to a spherical shell of thickness 4 = [1~3] km 
[7, 17]. Accordingly, we assume protons in the inner core, 
a sphere of radius i?i = i? — £s, to be in a normal (as op- 
posed to superconducting) phase. 

In a superconductor, components of the Maxwell stress 
tensor are enhanced by the ratio Hd/Bg, with respect to 



a normal conductor [47]. Here Hd ~ 10^^ G is the lower 
critical field [61] and Bg is the magnetic induction within 
the superconductor, which is always <^ Hd in a NS [62]. 
Several authors thus suggested [9, 48] that the deforma- 
tion caused by a typical magnetic field ^ 10^^ G would 
be ^ 10^ times larger than given by, e.g. Eq. (14). This 
would have very important implications for the scenario 
proposed here. 

For a phase transition in a pre-existing magnetic field, 
Ref. [46] argued that the magnetic induction in the su- 
perconducting outer core would be reduced with respect 
to that in the inner normal core, as required to maintain 
stress balance at the interface between the two regions. 
As a consequence, no significant change in the total NS 
ellipticity should be expected in this case. 

In the context of the present work, the toroidal mag- 
netic field is generated in the superconducting layer at 
the expenses of the r-mode energy. We thus expect 
that a larger deformation of the NS can be obtained 
only if a correspondingly larger amount of energy can 
be transferred from the modes to the toroidal field. Con- 
versely, for a fixed energy transfer rate, one would expect 
a weaker field to be produced in the superconductor, such 
that the total magnetic stress - thus the induced defor- 
mation - would be the same as for a normal conductor. 

A solution to the r-mode equations in an inhomoge- 
neous star would be required to solve this problem self- 
consistently [7], which is beyond our scope here. How- 
ever, we can gain some insight by considering the energy 
budget of the system. 

We adopt the expression of the magnetic damping rate 
in the normal and superconducting case, following [3] 
(cfr. Sec. HI). For the normal fiuid in the inner core, the 
rate of production of magnetic energy is (cfr. Eq. (31)) 






(''\t)^aUt)nit). (42) 



Here A and A are angular integrals of order unity defined 
in Ref. [3], $3.0 is the initial magnetic flux threading the 
inner core and B\ (t) is obtained by taking the volume 
average of the generated azimuthal field (cfr. Eq. (12)). 
The superconducting layer has a volume y^^*-^) — 
AnRlis and we can write 
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As the r-mode energy is tapped by the newly generated 
azimuthal field, the magnetic energy density and Maxwell 
stress grow both in the inner and outer core. We thus ex- 
pect that stress balance at the interface between the two 
regions will play a role in the subsequent development of 
r-modes. To obtain the rates of magnetic energy density 
production Ub in either region we divide Eqs. (42,43) by 



the corresponding volumes, y(^) and T/(^*^). Hence 
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having neglected numerical factors of order unity. If 
Maxwell stresses in the two zones need to meet an equilib- 
rium condition, the above ratio should be of order unity. 
Accordingly, a stable evolution of the internal magnetic 
field requires that asc/^N -- [Bf\t)/HciY/'^. This re- 
lation can be considered as a boundary condition on the 
amplitude at the surface separating the normal and the 
superconducting region. If this condition on the r-niode 
amplitudes is not satisfied, we expect that the magnetic 
instabilities will affect the growth of a in the two zones 
acting to restore the stress balance. From this, we con- 
clude that, if the condition of stress balance were to hold 
(at least in a time-average sense) , azimuthal fields of dif- 
ferent strengths would be generated in the two zones (see 
Eq. (44)), namely 
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This implies that the total NS deformation is expected 
to be of the same order of the deformation in a normal 
conductor. 

Note that superconducting protons may coexist with su- 
perfluid neutrons, in a fraction of the neutron star core. 
Pinning of neutron vortices to magnetic fiux tubes may 
affect in this case both the growth of the magnetic field 
in the superconducting shell and the spin frequency evo- 
lution [49, 50]. However, this complicated interplay and 
its ultimate implications are still subject of much debate 
[51-53] and the problem is far from settled. The coexis- 
tence of the two mutually interacting superfluid and su- 
perconducting particles of different species is still contro- 
versial, and notable observational limits on it have been 
proposed [54]. A detailed analysis of these effects and 
their inclusion in our numerical calculations are beyond 
the scope of this paper. 



Discussion 

The analysis of the previous scenarios allows us to draw 
some general conclusion about the internal magnetic 
fields and the maximum spin frequencies of the fastest 
accreting NSs. These conclusions are quite independent 
of the particular evolutionary path of the recycled stars. 
In the cases where i?0,ort ^ 10^'' G and es ~ Cmin, 
the GW emission due to the magnetic deformation does 
not limit the spin frequency of the accreting stars that 
can re-enter the r-mode instability window and gener- 
ate new azimuthal field. The new toroidal component 
can reach strengths B^ > 10^^ G causing a new mag- 
netic instability. On the other hand, in the cases where 
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FIG. 6: Limits on the maximuni spin frequencies of accreting 
neutron stars with an ellipticity |es| ~ 10~* due to the mag- 
netic deformation (solid line) . We show also the fastest accret- 
ing millisecond pulsars (blue points) and the burst oscillation 
sources (red points) [55]. The stars can enter horizontally in 
the gray area when the mass accretion rate decreases down to 
values M < lO~^M0yr~^ (see Discussion in Sec. VI). 



-B^.ort ^ 10^^ G the magnetic field rapidly decays by am- 
bipolar diffusion to strengths again of order 10^"' G. It is 
clear that an equilibrium configuration requires a value 
of the internal magnetic field -B^^ort ~ 10^'* G with an 
jesl ~ 10~^. This internal configuration limits the spin 
frequencies at values fmax ^ 630 MJg (e^/lO"^)"^/^ Hz 
preventing the star from further increasing their mag- 
netic field by re-entering the r-mode instability region. 
In Fig. 1 and Fig. 6 we show the limits on the maximum 
spin frequencies of accreting neutron stars as a function 
of temperature and of mass accretion rate respectively. 
In Fig. 6 we plot also the observational data of accreting 
millisecond pulsars and of burst oscillation sources. The 
fastest spin frequencies {v > 650 Hz) can be reached only 
for mass accretion rates M > lO~^M0yr~-'^ with a spin 
up time-scale in the range « [10^-10^] yrs. Note that 
during this period of time, the internal magnetic field 
B^^ort is not subject to significant decay. 

Taking into account the maximum values of M for the 
NSs in LMXBs (w 10~^ Mq yr~^) we can estimate a max- 
imum spin frequency Vmax ^ 800 Hz which is in agree- 
ment with the estimate lymax ^ 730 Hz of Rcfs. [10, 11]. 

Note that the star can move horizontally in Fig. 6 en- 
tering the forbidden region (gray area) when the mass ac- 
cretion rate decreases down to values M < lO^^Moyr"^. 
An accreting spinning up neutron star is expected to 
move upwards in the white area in Fig. 6. When its 
spin frequency reaches the limiting value indicated by 
the solid line the star cannot spin up further and its po- 
sition remains close to the limiting line. If subsequently 
the mass accretion diminishes, the star drift horizontally 
to the left Fig. 6. Therefore, in this interpretation, stars 



10 



lying in the gray area of Fig. 6 are presently accreting at 
significant lower rate than in the past. 



VII. DETECTABILITY OF THE 
GRAVITATIONAL- WAVES EMITTED 

In this section we discuss the detectability of the grav- 
itational radiation emitted by recycled millisecond NSs 
due their magnetically-induced distortion. We consider 
a typical distance d = 10 kpc of accreting NSs in our 
Galaxy and we calculate the average gravitational wave 
amplitude h. 
The instantaneous signal strain can be expressed as [21]: 
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where B^^is = ^^/(lO" G), dxn = d/(10kpc) and z^soo = 

://(500Hz). 

The spin frequency of accreting neutron stars changes 

significantly over a (long) timescale 
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where Km = {G/MR'^Y^'^ /I. Hence we can integrate for 
long periods Tabs- The minimal detectable signal ampli- 
tude ho is [55] 



ho « 11.4 



t^n 



(48) 



where Sn is the power spectral density of the detector 
noise. 

In Fig. 7 we show h as a, function of the spin frequencies 
i/ of the star, while in Fig. 8 we compare the predicted 
and detectable amplitudes as a function of the frequency 
f = 2v oi the signal. 

We show also the sensitivity curve of LIGO and Virgo 
and of the next generation of detectors (Advanced LIGO, 
Advanced Virgo and Einstein Telescope). 



VIII. CONCLUSIONS 

We have shown that accreting millisecond NSs can be 
deformed significantly by the very large magnetic fields 
generated by r-modes. A secular instability takes place 
when the magnetic distortion dominates over the max- 
imal oblateness that the NS crust can sustain (i?^ > 
10^'^ G). As a consequence of this instability, the angular 
momentum J* and the magnetic axis become orthogonal 
on a timescale Tv given by Eq. (17) and the star begins 
to emit GWs. 

We have shown that the GWs emission due to the 
magnetic deformation can limit the spin frequencies of 
accreting NSs. In particular we obtain a maximum spin 



frequency Vmax ^ 800 Hz. Our results are quite inde- 
pendent of the exact evolution of the internal magnetic 
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FIG. 7: Averaged gravitational wave amplitude h vs. the 
spin frequency of the star v. The gray area is obtained as- 
suming that the accreting neutron star is located at a distance 
d = 10 kpc and that the internal magnetic field is in the range 
B^ ~ [0.2-1] X 10^'* G and fc^ — 1. For comparison we show 
also the maximum amplitude of the GWs emitted by the star 
due to r-modes (dotted-dashed line) after the rotation of the 
internal magnetic field or the development of the Tayler in- 
stability. We assume that a ~ 10~* (see Figs. 4 and 5). 
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FIG. 8: We compare the predicted amplitude (gray area) and 
the detectable amplitude for LIGO (2010) [56], Virgo [57], Ad- 
vanced LIGO [58] , Advanced Virgo [59] and the Einstein Tele- 
scope [60] as a function of the frequency of the gravitational 
waves. Here we consider an integration time T^bs ~ 1 yr. As 
in Fig. 7, d = 10 kpc, B^ « [0.2-1] x 10^" G and k, = 1. 



field, that is difficult to estimate and depends on several 
unknown parameters, e.g. the growth timescale of the 
Tayler instability and the parameter n that is related to 
the flip time of the generated toroidal field. 

In the end we analyzed the GWs emission due to the 
magnetic distortion. Our results suggest that with an in- 
tegration time of 1 year, the next generation of detectors 
(e.g. the Einstein Telescope) should be able to detect 
GW signals by accreting millisecond NSs located in our 
Galaxy. 
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This is the minimal magnetic field strength required to 
force a non-zero magnetic fiux through a type-II super- 
conductor. 

Strictly speaking, the NS core should be in the Meissner 
state with total flux expulsion. However, the large elec- 
trical conductivity of the core material forces it into a 
metastable type-II state, cfr. [45] 



